On Electron Transport in ZrBi 2 , ZrB 2 and MgB 2 
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We report on measurements of the temperature dependence of resistivity, p(T), for single crystal 
samples of ZrBi2, ZrB2 and polycrystalline samples of MgB2. It is shown that cluster compound 
ZrBi2 behaves like a simple metal in the normal state, with a typical Bloch - Griineisen p(T) 
dependence. However, the resistive Debye temperature, Tr — 300 K, is three times smaller than To 
obtained from specific heat data. We observe the T 2 term in p(T) of these borides, which could be 
interpreted as an indication of strong electron-electron interaction. Although the p(T) dependence of 
ZrBi2 reveals a sharp superconductive transition at T c = 6.0 K, no superconductivity was observed 
for single crystal samples of ZrB2 down to 1.3 K. 

PACS numbers: 74.70.Ad, 74.60.Ec, 72.15.Gd 



It is known that boron has a tendency to form cluster 
compounds. In particular there are octahedral Bg clus- 
ters in MeBg, icosahedral B12 clusters in /?-rhombohedral 
boron, and cubo-octahedral B12 clusters in MeBi2. So 
far, several superconducting cubic hexa - MeB@ and do- 
decaborides - McBi 2 have been discovered [lj (Me=Sc, 
Y, Zr, La, Lu, Th). Many other cluster borides (Me= Ce, 
Pr, Nd, Eu, Gd, Tb, Dy, Ho, Er, Tin) were found to be 
ferromagnetic or antiferromagnetic [j], |2| ■ Even though 
the superconductivity in ZrBi2 was discovered a long 
time ago (T c = 6 K) [l|, there has been little effort de- 
voted to the study of electron transport and basic super- 
conductive properties of dodecaborides. Only recentl 
the electron transport of solid solutions Zri_ a; Sc a ;Bi2 
as well as the band structure calculations of ZrBi2 4] has 
been reported. Understanding the properties of the clus- 
ter borides as well as the superconductivity mechanism 
in these compounds is very important. 

Recently, we reported superconductivity at 5.5 K in 
the polycrystalline samples of ZrB2 [5|. This was not 
confirmed in later studies |(j- It was recently suggested 
y| that this observation could be associated with nonsto- 
ichiometry in the zirconium sub-lattice. In this letter we 
address this problem. We present the results from mea- 
surement of the temperature dependencies of resistivity, 
p(T), for single crystals of ZrBi2 and ZrB2. Compara- 
tive data from polycrystalline samples of MgB2 are also 
presented. The superconducting properties of ZrBi2 will 
be published elsewhere. 

Under ambient conditions, dodecaboride ZrBi2 crys- 
tallizes in the fee structure of the UB12 type (space 
group Fm3m), a = 0.7408 nm 8]. In this structure, 
the Zr atoms are located at interstitial openings in the 
close-packed B12 clusters 3]. In contrast, ZrB2 shows a 
phase consisting of two-dimensional graphite-like mono- 
layers of boron atoms with a honeycomb lattice structure, 



intercalated with Zr monolayers (with lattice parameters 
a = 0.30815 nm and c = 0.35191 nm @). 

The ZrB2 powder was produced by the boron carbide 
reduction of Zr02- The ZrBi2 single crystals were ob- 
tained from a mixture of a certain amount of ZrB2 and 
an excess of boron (50 — 95%). The resulting mate- 
rials were subjected to a crucible-free RF-heated zone- 
induction melting process in an argon atmosphere. The 
obtained single crystal ingots of ZrBi2 and ZrB2 have 
a typical diameter of about 5 — 6 mm and a length of 
40 mm. A metallographic investigation detected that 
the ZrB2 crystal is surrounded by a polycrystalline rim 
about 0.5 mm thick. The measured specific density of 
the ZrBi2 rod is 3.60 g/cm 3 , in good agreement with the 
theoretical density. The X-ray diffraction measurements 
confirmed that both ingots are single crystal. We found 
the cell parameters of ZrBi 2 , a = 0.74072 ± 0.00005 nm, 
to be very close to the published values p|. 

Polycrystalline MgB2 and CaMgB2 samples were sin- 
tered from metallic Mg or a mixture of Ca,Mg powders 
and boron pellets using a similar technique as outlined in 
our earlier work [5j . This technique is based on the reac- 
tive liquid Mg,Ca infiltration of boron. X-ray diffraction 
patterns and optical investigation show large grains of 
single MgB2 phase, with much smaller grains of semi- 
conducting CaBg phase visible in-between. Density of 
MgB2 grains was rather high, 2.4 g/cm 3 , while the sam- 
ples prepared from Mg infiltration had smaller density 
of 2.2 g/cm 3 . Only MgB2 samples cut from large grains 
were studied. These samples will be denoted as CaMgB2. 

We used a spark erosion method to cut the samples 
into a parallelepiped with dimensions of about 0.5 x 
0.5 x 8 mm 3 . Single crystal samples were oriented along 
<100> for ZrBi2, and in hexagonal [0001] and basal 
[1100] directions for ZrB2, respectively. The orientation 
process was performed using an X-ray Laue camera. The 
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FIG. 1: Temperature dependence of the resistivity, p(T), 
of ZrBi2 single crystal (open circles), MgB2 (squares) and 
CaMgB2 (triangles) samples. The solid lines represent BG 
fits to the experimental data by Eq. 



samples were lapped by diamond paste and subsequently 
etched: ZrBi2 in hot nitrogen acid, ZrB2 in mixture of 
H2O2/HNO3/HF, and MgB 2 in 2%HC1 plus water-free 
ethanol. 

A standard four-probe ac (9 Hz) method was used for 
resistance measurements. We used Epotek H20E silver 
epoxy for electrical contacts. The samples were mounted 
in a a temperature variable liquid helium cryostat. Tem- 
perature was measured with platinum (PT-I03) and car- 
bon glass (CGR-I-500) sensors. The critical temperature 
measured by RF susceptibility [f| and p(T) was found to 
be T c0 = 5.97 K for ZrBi 2 samples and 39 K for MgB 2 
polycrystallinc samples, respectively. 

We display the temperature dependence of the resis- 
tivity for ZrBi2, MgB 2 and CaMgB 2 in Fig. [JJ and that 
of ZrB2 in Fig. [3 To emphasize the variation of p(T) in 
a superconductive state, we plot these data in the inset 
of Fig. ^ The samples demonstrate a remarkably narrow 
superconducting transition with AT = 0.04 K for ZrBi2 
and with AT = 0.7 K for both MgB 2 samples. Such a 
transition is a characteristic of good quality samples. 

As we can see from Fig. |5J no superconductivity was 
observed in ZrB 2 down to 1.3 K, while a pronounced 
slope change in p(T) is observed around 7 K. One can 
explain such behavior in the following way. In ZrB2 the 
Fermi level is located in the pseudo gap. The presence 
of Zr defects in Zro.7sB2 leads to the appearance of a 
very intense peak in the density of states in the vicinity 
of the pseudo gap and subsequent superconductivity Q • 
We strongly believe that the observation of y| was due 
to nonstoichiometry of our samples. Superconductivity 
in nonstoicheometric samples is very common in other 
borides: MoB 2 . 5 , NbB 2 . 5 , Mo 2 B, W 2 B, BeB 2 . 7 5 SIS- 



FIG. 2: Temperature dependence of p(T) of ZrB2 single 
crystal samples in basal plane (circles) and in c direction 
(squares). 



It is worth noting that ZrBi2 is mostly boron, and 
one could speculate that its resistivity should be rather 
high. In contrast we observe that the room temperature 
resistivity of ZrBi2 is almost the same as for MgB2 and 
ZrB 2 samples. The p(T) is linear above 90 K with the 
slope of p(T) more pronounced than in MgB2 or ZrB2. 
The residual resistivity ratio RRR of 9.3 for ZrBi2 as 
well as RRR « 10 for MgB2 and ZrB2 samples suggests 
that the samples are in the clean limit. One can predict 
a nearly isotropic resistivity for fee ZrBi2, which can be 
described by the Bloch-Grunciscn (BG) expression of the 



electron-phonon e-p scattering rate 
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Here, p(0) is the residual resistivity, p\ = dp(T)/dT 
is a slope of p(T) at high T (T > T R ), t = T/T Rl T R is 
the resistive Debye temperature and J$(l/t) is the Debye 
integral. As we can see from Fig.Q ah data for ZrBi2 fall 
very close to the theoretical BG function (solid line). To 
emphasize the variation of p{T) at low T we plot these 
data as p(T) — p(0) versus t 5 J 5 (l/t) in Fig.|3]on a log-log 
scale. The BG formula predicts a linear dependence of 
log\p(T) — p(0)] versus log[t 5 J 5 (l/t)] with the slope equal 
to unity. We use T R as a fitting parameter to achieve 
agreement at the high temperatures. For comparison, we 
also present our p(T) data of ZrB2 and MgB2 calculated 
in a clean case of the two band model [12| • 

It is clear from Fig. [3] that above 25 K the BG model 
describes the p(T) dependence of ZrBi2 fairly well. It is 
remarkable that this description works well with constant 
T R = 300 K. At the same time, To calculated from 
specific heat data [13J is three times higher. Furthermore 
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FIG. 3: The p(T) - p(0) vs. reduced Debye integral f 5 J 5 (l/i) 
for ZrBi2 (open circles), ZrE>2 in basal plane (crosses) and 
CaMgE<2 (squares). The dashed line is p(T) of MgB2 calcu- 
lated in the two band model |l2| . 



Td increases from 800 K to 1200 K as temperature varies 
from T c up to room temperature. In order to shed light 
on this discrepancy, we used a model applied to LaBt; of 
Ref. [lj| . We can treat the boron sub-lattice as a Debye 
solid with Tji and the Zr ions as independent Einstein 
oscillators with characteristic temperature Tg. The effect 
of the Einstein mode on the resistivity of a metallic solid 
is discussed in Ref. 15]: 



Pe(T) 
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Here N is the number of oscillators per unit volume, 
if is a constant that depends on the electron density of 
the metal, M is the atomic mass. We fit the data by 
summing Eq. (JTJ and Eq. @, and living KN/M, pi, T R 
as free parameters. Although the model calculations per- 
fectly match the data (see solid line in Fig.[3|), the Te we 
are getting is unreasonably small (Te = 50 K), and the 
difference between Tr and specific heat To becomes even 
worse, T R = 270 K. We believe that this in-consistency 
of Tr and To can be explained by limitation of Tr by 
a cut-off phonon wave vector q — k R T /hs. The latter 
is limited by the Fermi surface (FS) diameter 2Uf [llj 
rather than the highest phonon frequency in the phonon 
spectrum. 

According to band structure calculations [4J, the FS 
of ZrBi2 consists of an open sheet along TL direction at 
point T with krx — 0.47 A -1 , a quasi spherical sheet at 
point X (kxr = 0.37 A^ 1 ) and a small sheet at point 
K (k Kr = 0.14 A" 1 ). We suggest that T R is limited 
by the small FS sheet. Unfortunately the experimen- 
tal FS model and the sound velocity are not yet known. 



Therefore we can not corroborate this suggestion by ex- 
perimental FS. 

As we can see from Fig. [31 the p(T) of ZrB2 and MgB2 
samples deviates from the BG model even more dramat- 
ically. Putti et al. 17] modified the BG equation intro- 
ducing variable power n for the t n J n (l/t) term in Eq. QJ. 
The best fit to the data was obtained with n = 3 which 
in fact ignores a small angle e-p scattering. Recently So- 
logubenko et al. [lg] reported a cubic T dependence in 
the a,b plane resistivity below 130 K in the single crys- 
tals of MgB2. This was attributed to the interband e-p 
scattering in transition metals. 

However, we believe there are strong objections to this 
modified BG model: (i) a cubic p(T) dependence is a the- 
oretical model for large angle e-p scattering and no evi- 
dence of it was observed in transition and non-transition 
metals; (ii) the numerous studies of the p(T) dependence 
in transition metals have been found to be consistent 
with a sum of electron-electron e-e, T 2 , and e-p, T 5 , con- 
tributions to the low T resistivity, which may easily be 
confused with a T 3 law [UdES]; ( ui ) the interband 
a — tt e-p scattering plays no role in normal transport in 
the two band model for MgB2 [lg . 

In order to solve these problems, we added e-e scat- 
tering T 2 term in Eq. QJ 19, 20] as a possible scenario. 
Indeed, keeping in mind that the BG term is proportional 
to T 5 at T < O.ITr, p(T) dependence may be presented 



in a simple way 



[p(T) - p(0)]/T 2 = a + (3T*. 



Here a and j3 — 497.6pi/T R are parameters of e-e and e- 
p scattering terms, respectively. Such a plot should yield 
a straight line with slope of (3 and its intercept with y- 
axis (T = 0) should equal to a. Further, to be consistent 
with BG law, the (3 parameter should lead to the same 
Te, as obtained from high T log-log fit in Fig. [31 and both 
coefficients must be independent of p(0). We determined 
p(0) from the intercept of linear p(T) vs T 2 dependence 
with the T = axis and plotted the [p(T) - p(0)]/T 2 vs. 
T 3 in Fig. 01 It is evident that the measured resistivity 
approaches a quadratic law at T < 25 K in ZrBi2, at 
T < 100 K in ZrB 2 , and at T < 150 K in both MgB 2 
samples. 

The regime of applicability of two term fit is limited 
to temperatures below O.lTfl. At larger T the e-p term 
increases more slowly than T 5 law and this is why the 
data are not consistent any more with the two terms 
equation. From the intercept with T — axis, we find 
very similar values of a for ZrBi2 and ZrB2 samples in 
the basal plane (a = 22 pflcmK~ 2 and 15 ptlcmK ~ 2 , re- 
spectively) while a is about five times larger for CaMgB2 
sample, 95 pVtcmK~ 2 . The slopes of j3 give p\ and Tr 
values largely consistent with high temperature log-log 
fits for the ZrBi2 and ZrB2 samples. 

However, low T results for (3 and p\ are far from con- 
sistent with high T data for both MgB2 and CaMgB2 
samples. Nevertheless, the magnitude of Tr = 900 K 
for MgB2 extracted from log-log fit above 150 K, is in 
excellent agreement with Tq = 920 K obtained from 
low-temperature specific heat measurements [2l|, and is 
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FIG. 4: Low temperature behavior of \p(T) — p(0)]/T 2 versus 
T 3 for: (a) ZrBi2 (circles), ZrB2 in basal plane (squares), 
ZrB2 along c (triangles) and (b) MgB2 (circles) and CaMgB2 
(squares) samples. 



considerably lower then the reported data based on T 3 
dependence of p(T)(Tr = 1050-1226 K, where T 2 term 
was ignored. [g, [171 llSjL A similar fit for theoretical 
curve is even more consistent with Tr = 900 K, however 
we have to mention that violation of Matthiessen's rule 
in MgB2 may mask the intrinsic p(T) dependence [l2fl. 

In general, there are many scattering processes respon- 
sible for the T 2 term in p{T) of metals: (i) size, surface, 
dislocation and impurity scattering induced deviations 
from Matthiessen's rule (see references in [22|); (ii) e- 
p scattering for small cylindrical FS sheets relative to 
the phonon wave vect or ll6J : (iii) inelastic electron im- 
purity scattering (e-i) J23J; (iv) t he q uantum interference 
between e-i and e-p scattering J24(; (v) e-e scattering 

urn. 

We can estimate some of these effects. We use Drude 
law to obtain the residual electron mean free path I = 
Attvf/p^p- Using a Fermi velocity of v„ = 3.2 • 10 7 cm/s 
and a plasma frequency cjp = 5.16 • 10 15 s _1 for MgB 2 
CT-band [1^. we obtain I w 100 nm. This implies that 
size effects are negligible for both MgB2 samples and Zr 



borides. In agreement with ZrB 2 data (see Fig. Bt the a 
is proportional to p(0) for inelastic e-i scattering |23ll24|| . 
However, this term is 1.5 times lower for CaMgB2 relative 
to MgB2, which has the same p(0). 

We can try to estimate contribution from the small 
FS sheets to a. The T 2 term was observed in p(T) and 
electron scattering rates of Bi and Sb, which was at- 
tributed to a missing of one q component for e-p scat- 
tering on small cylindrical FS sheets [Ig. The FS of 
MgB2 is composed of two warped open cylinders run- 
nin g al ong the c axis, which arise from a boron orbitals 
[12I l2q . The FS of ZrB 2 consist of nearly ellipsoidal sur- 
faces joined together at the corners J2a, |22fl , which may 
also be responsible for the T 2 term in p(T). We can use 
the sound velocity s = 1.1 • 10 6 cm/s and 8 • 10 5 cm/s 
for MgB2 and ZrB2, respectively [28|, |2!j, to estimate 
the lowest temperature, T min = fikps/ks, when the 
phonon wave vector q matches a neck of smaller a tube 
in MgB 2 (k a = 0.129 A" 1 [25|) or a diameter of the ellip- 
soidal sheets in ZrB 2 (k F = 0.095 A" 1 ) |2fj). We obtain 
T m in = 95 K and 60 K, respectively. Thus we conclude 
that q < kp at T < 100 K in both diborides, which im- 
plies that the contribution of the 2D FS sheets to a is 
negligible. 

In general only umklapp e-e scattering contributes to 
p{T), whereas the normal collisions are significant in com- 
pensated metals and in thermal resistivity p0| . Borides 
have rather high T^> which depresses the e-p scattering, 
so that the e-e SR term is easier to observe. Notice how- 
ever that the a value for MgB2 is five times larger than 
corresponding values in ZrBi2 and ZrB2. The latter val- 
ues are in turn five times larger than in transition met- 
als {oLMn = 2.5 pVtcm/K 2 and aw — 1.5 — 4 pilcm/K 2 
[T9I l20lp. Therefore, additional experiments must be 
performed for more pure samples before final conclusion 
about the origin of the T 2 term in borides can be drawn. 

In conclusion, we present a study of the p(T) of single 
crystals of ZrBi2, ZrB2 and polycrystalline samples of 
MgB2. Large differences between resistive and specific 
heat Debye temperatures have been observed for ZrB i2 . 
The results provide evidence of a T 2 term for all these 
borides at low T, whose origin is not yet understood. 
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